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ABSTRACT 

This  report  summarizes  analytical  studies  conducted  during 
the  past  year  to  determine  the  interaction  of  sonar  domes  and 
transducers  during  transmission.  A  general  description  of  the 
mathematical  model  used  in  the  analysis  is  given  together  with 
some  improvements  which  have  been  incorporated  in  the  model  dur¬ 
ing  the  past  year.  The  results  of  a  study  of  the  dependence  of 
certain  performance  criteria  such  as  directivity  index  on  the 
parameters  of  the  dome -transducer  model  is  reviewed.  In  addition, 
a  general  description  of  developments  related  to  the  treatment 
of  viscoelastic  dome  materials  and  nonconcentric  dome  shapes 
is  given. 

The  conclusions  drawn  in  the  various  technical  memoranda 
submitted  during  the  year  are  restated  in  this  report.  Briefly 
these  are:  (a)  the  dome  tends  to  offset  the  gains  in  side  lobe 
reduction  obtained  by  shading  the  transducer  elements,  and  (b) 
the  most  important  property  of  the  dome  model  in  determining  its 
effects  on  various  transducer  performance  criteria  is  its  surface 
density  (product  of  mass  density  and  thickness).  These  conclusions 
are  supported  by  numerical  results  presented  in  the  report. 

Finally,  some  recommenda cions  for  the  continuance  of  studies 
of  dome -transducer  interactions  are  given.  These  include  an 
emphasis  on  the  development  of  a  dome -transducer  model  of  more 
general  configuration  and  the  desirability  of  an  experimental 
program  to  support  the  analytical  studies. 
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I.  INTRODUCTION 

This  report  is  a  summary  of  the  developments  achieved 

during  the  past  year  in  a  program  of  studies  to  analyze  the 

interactions  of  sonai  domes  and  transducers  during  transmission. 

The  program  of  analysis  has  been  based  on  an  analytical  model 

consisting  of  a  circular  cylinder  (transducer)  and  concentric 

shell  (dome) .  The  basic  model  was  developed  under  previous  studies 

.  x 

and  reported  in  the  summary  report  for  C.Y.  1964. 

At  the  beginning  of  the  past  year's  study,  it  was  decided 
that  a  thorough  understanding  of  the  existing  model  would  lead 
to  some  elementary  guidelines  for  improved  dome  design.  Accord¬ 
ingly,  some  improvements  were  made  in  the  existing  model,  speci¬ 
fically,  an  improved  description  of  the  velocity  distribution  on 
the  transducer  face;  and  an  evaluation  of  the  dependence  of  some 
well-known  performance  criteria  on  the  parameters  of  the  dome- 
transducer  model  was  conducted. 

As  a  result  of  these  studies  it  is  now  possible  to  make 
qualitative  statements  about  the  relative  importance  of  the 
system  parameters  in  determining  the  effect  of  the  sonar  dome  on 
system  performance.  For  example,  it  is  known  that  the  surface 
density  of  the  dome,  i„e„,  the  product  of  the  dome  thickness  and 
dome  material  density,  is  the  most  important  parameter  in  deter¬ 
mining  the  acoustic  properties  of  the  dome.  The  stiffness 
properties  (elastic  moduli)  of  the  dome  material  are  relatively 
unimportant  in  determining  acoustic  properties.  It  has  been 
shown  that  the  spacing  of  the  dome  relative  to  the  transducer 
can  be  a  significant  factor  in  determining  the  radiation  impedance 
at  the  transducer  face  and  the  maximum  energy  transferred  to  the 
medium  by  the  transducer. 

These  studies  have  given  rough  indications  for  improve¬ 
ments  in  dome  design.  However,  analytical  models  of  the  dome- 
transducer  must  be  considered  as  still  in  the  primitive  stage. 
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Analytical  and  related  experimental  studies  must  be  pursued  in  a 
continuing  exploratory  development  program  in  order  to  achieve 
further  improvements  in  sonar  dome  design. 
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II.  BACKGROUND 

The  proper  design  of  a  sonar  dome  must  include  several 
aspects  including  hydrodynamic,  strength,  durability,  and.  of 
course,  acoustic  properties.  While  the  dome  can  be  studied 
independently  in  terns  of  those  requirements,  the  final  design 
must  be  a  marriage  of  the  results  of  the  various  studies  and 
must,  therefore,  represent  a  trade-off  among  each  of  the  require¬ 
ments.  The  studies  conducted  at  TRACOR  have  been  directed 
principally  toward  the  acoustical  aspect  of  dome  design.  The 
broad  purpose  of  these  studies  is  to  define  guidelines  for 
improved  sonar  performance  which  can  be  integrated  into  a  total 
design  concept  for  domes  including  hydrodynamic,  strength,  and 
durability  requirements . 

In  order  to  develop  an  understanding  of  the  acoustical 
implications  of  a  dome,  it  is  helpful  to  consider  a  simple  two 
dimensional  cylindrical  transducer.  The  sound  pressure  p  of 
this  transducer  can  be  represented  as 

H^(kr)  cos  n0  e  ,  (1) 

where  H^(kr)  is  the  first  type  of  Hankel  function  of  argument  lcr 

and  order  n,  k  =  ~  with  iu  the  angular  frequency,  c  is  the  speed 

of  sound  in  the  medium,  r  and  0  are  polar  coordinates ,  a  is  the 

transducer  radius ,  and  the  prime  indicates  a  derivative  with 

respect  to  the  argument  of  the  Hankel  function.  Steady  state 

♦  . 

time  dependence  of  the  form  e”  w  is  assumed.  The  6  are  Fourier 
coefficients  determined  by  a  specification  of  the  normal  particle 
velocity  on  the  transducer  face.  Generally,  the  velocity  speci¬ 
fication  is  such  that  the  radiating  portion  of  the  transducer  is 
divided  into  segments  having  plane  wave  phasing  or  time  delay 
and  an  arbitrary  shading,  for  example,  cosine  shading. 


=  2 


n 


n==0  n 


(ka) 


3 


TRACOR,  INC. 


6500  TRACOR  LANE,  AUSTIN.  TEXAS  78721 


Equation  1  can  be  used  to  compute  nearfield  and  farfield 

sound  pressure  distributions,  and  subsequently  determine  a  two- 

dimensional  directivity  index,  major  lobe  width,  and  side  lobe 

level.  The  sound  pressure  at  any  point  in  space  is  determined 

approximately  by  a  finite  number  of  terms  or  "modes,"  the  number 

depending  on  ka,  kr  (where  r  is  the  point  of  interest),  and  the 

form  of  6  ■. 

n 

If  the  same  cylindrical  radiator  is  surrounded  by  a  con¬ 
centric  elastic  shell,  the  expression  describing  the  pressure, 

2 

p^ ,  outside  the  shell  (See  Fig.  1.)  is 


6 

n 


(ka)  +  ZncinCjn(kb)  tT(ka)  -  N^kb)  J^ka)] 


Hb(kr)  cos  n9  e  ia,t  ,  (2) 

where  b  is  the  mean  radius  of  the  shell,  J  and  Np  are  the  Bessel 
and  Neumann  functions  respectively,  =  (iu)nb/2ocZ)Hn(kb)  ,  zn  is 
referred  to  as  the  modal  impedance  and  the  remaining  symbols  are 
as  previously  defined.  The  shell  now  can  be  considered  a  new 
radiator  which  produces  the  sound  field  given  by  Eq.  2.  If  .  ie 
sound  pressures  p  and  P2>  Eqs .  1  and  2, are  compared  at  some  point 
in  space  outside  the  shell  (fixed  r  and  0),  they  will  not,  in 
general,  be  the  same.  Hence,  the  dome  affects  the  sound  field 
produced  by  the  transducer.  In  order  that  the  dome  not  affect 
the  transmitted  beam  it  is  necessary  that  Zn  be  zero  for  all 
values  of  n  ([J^(kb)  N^(ka)  -  Nn'(kb)  J^(ka)]  in  general  is  non¬ 
zero  for  a>0 ,  b>0) .  In  general  Zn  is  nonzero  for  all  except 
one  or  two  values  of  n,  i.e.,  one  or  two  terms  in  Eq.  2.  Hence, 
for  the  same  velocity  distribution  on  the  transducer,  the  pressure 
field  at  some  distance  from  the  transducer  will  be  different  from 
that  of  a  bare  transducer  if  the  transducer  is  surrounded  by  a 
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concentric  shell.  Computations  of  farfield  directivity  patterns 
indicate  that  the  presence  of  the  shell  results  in  an  increase 
in  the  level  of  the  side  lobes  and  a  reduction  of  directivity 
index  compared  to  a  bare  transducer. 

The  shell  also  affects  the  sound  pressure  distribution 
near  the  transducer.  The  general  form  of  the  sound  pressure  p^ 
between  the  shell  and  transducer  is 

cc 

P1  =  2  ®-An  Hn^kr)  +  BnHn(kr)3  C0S  nG  »  O) 

n=0 


where  A  and  B  are  constants  that  depend  on  6  and  Z  ,  and 
H^(kr)  is  the  Hankel  function  of  the  second  kind.  The  presence 
of  both  types  of  Kankel  functions  indicates  that  there  are  acous¬ 
tic  waves  traveling  both  outward  from  the  transducer  and  inward 
toward  the  transducer.  The  solution  form  for  a  bare  transducer, 

Eq.  1.  contains  only  outward  traveling  waves.  The  inward  travel¬ 
ing  waves  represent  reflections  from  the  dome  which  arise  due  to 
the  impedance  discontinuity  of  the  dome.  If  Zn  were  zero  for 
each  value  of  n,  then  Eq.  3  would  reduce  to  Eq.  1  with  An  = 

6 

-r-r~ -  and  B  =0.  Again,  however,  Z  can  be  zero  for  only  one 

H;  (ka)  n  n 

n 

or  two  values  of  n,  and  the  presence  of  the  shell  will  affect  the 
sound  pressure  distributions  near  the  transducer  compared  to  a 
transducer  with  no  shell.  Computations  of  nearfields  for  trans¬ 
ducers  with  and  without  shells  indicate  that  the  shell  significantly 
alters  the  radiation  loading  on  the  transducer.  In  addition,  the 
nearfield  sound  pressure  distribution  for  a  transducer  with  a 
shell  is  highly  nonuniform  compared  to  the  nearfield  of  a  bare 
transducer. 

The  present  model  assumes  the  dome  to  be  a  homogeneous 
body  with  no  stiffening  members.  Intuitively  one  feels  that  the 
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addition  of  stiffening  structures  to  domes  will  tend  to  degrade 
their  performance.  The  effect  of  stiffening  can  be  visualized 
by  considering  that  the  stiffening  members  tend  to  divide  the 
dome  up  into  a  series  of  plates  whose  edges  are  constrained.  The 
edge  constraints  impose  a  somewhat  different  vibration  response 
of  the  dome  compared  to  an  unstiffened  shell.  As  a  result,  it 
is  expected  that  the  stiffening  would  further  distort  the  farfield 
pressure  distribution.  A  comparison  of  some  acoustic  characteris¬ 
tics  of  stiffened  and  unstiffened  plates  has  been  carried  out 
3 

elsewhere,''  and  the  indication  is  that  the  effect  of  stiffening 
structure  is  significant,  particularly  in  terms  of  increased 
transmission  loss.  Unfortunately,  the  adaptation  of  analytical 
methods  used  in  stiffened  plate  studies  to  the  dome -transducer 
model  is  not  feasible  at  this  time. 

The  studies  conducted  to  date  have  given  some  insight 
into  the  problem  of  improving  the  acoustic  performance  of  domes. 
While  the  dome  impedance  cannot  be  made  zero  for  all  modes  of 
vibration,  it  is  possible  to  reduce  the  modal  impedance  for  all 
modes  and  thereby  improve  the  acoustic  properties  of  the  dome 
(at  least  for  unstiffened  domes) .  The  methods  for  obtaining 
this  improvement  are  described  in  the  following  sections. 
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III.  DISCUSSION  OF  ANALYTICAL  DEVELOPMENTS 

This  section  of  the  report  contains  descriptions  of  the 

analytical  developments  in  the  dome -transducer  model  during  the 

past  year.  Some  of  these  developments  have  been  reported  in 

2  4 

separate  technical  memoranda,  ’  and  the  reader  is  referred  to 
these  for  details  of  the  analyses. 

A.  Improvement  of  the  Transducer  Boundary  Specifications 

As  described  in  Section  II,  the  effects  of  the  sonar 
dome  on  transmitted  beams  have  been  computed  with  a  mathematical 
model  consisting  of  a  cylinder  (transducer)  and  concentric  shell 
(dome) .  To  form  a  typical  KDT  beam,  a  portion  of  the  transducer 
elements  must  be  driven  in  a  prescribed  manner  —  usually  the 
elements  are  phased  or  time-delayed  to  straight  line  and  have 
amplitude  shading.  In  the  first  analytical  model  of  the  dome- 
transducer,  it  was  not  possible  to  describe  the  transducer  element 
velocities  with  discrete  amplitudes  and  phases.  It  was  necessary 
to  describe  the  velocity  amplitude  and  phase  distributions  on  the 
transducer  face  as  continuous  functions.  The  first  improvement 
in  the  dome -transducer  model  was  the  development  of  the  procedures 
necessary  to  specify  discrete  values  of  velocity  amplitude  and 
phase  for  each  of  the  elements  used  in  forming  a  transmit  beam. 

(The  details  of  this  improvement  are  given  in  Ref.  4.). 

This  improvement  results  in  several  advantages.  First, 
the  model  is  a  more  realistic  representation  of  the  physical 
device  and  one  expects  the  results  of  the  model  to  predict  more 
accurately  the  phenomenon  being  studied.  Second,  the  effect  of 
the  dome  on  transmitted  beams  for  various  types  of  element  shad¬ 
ing  can  be  predicted.  Third,  the  model  can  be  used  to  predict 
the  effect  of  inoperative  elements  on  farfield  directivity 
patterns  (by  assigning  zero  velocity  amplitude  to  the  inoperative 
elements)  and  thereby  establish  tolerance  limits  on  system 
performance . 
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Comparisons  of  results  for  the  older  dome -transducer 
model  and  the  newer  model  with  improved  transducer  boundary 
specification  were  made  in  Ref.  4.  Significant  differences  in 
both  farfield  and  nearfield  sound  pressure  distributions  were 
obtained  using  the  two  models.  (See  Figs.  2  and  3.)  The  newer 
model  also  was  used  to  determine  the  farfield  sound  pressure 
distributions  for  several  types  of  element  shadings.  Results 
for  cosine  shading  and  Doiph-  Tchebycheff  shading  are  shown  in 
Figs.  4  and  5.  These  studies  indicated  that  increased  side  lobe 
level,  the  primary  effect  of  the  dome  on  the  farfield  pattern, 
was  more  pronounced  for  cases  of  severe  shading  such  as  Dolph- 
Tchebycheff . 

B .  Parametric  Analysis  of  Dome -Transducer  Interactions 

More  recently,  the  improved  dome -transducer  model  has 

been  used  to  compute  dome  effects  for  a  systematic  variation  of 

2 

the  model  parameters.  In  particular,  the  dependence  of  such 
performance  criteria  as  directivity  index,  side  lobe  level,  major 
lobe  width,  and  transmission  coefficient  on  dome -transducer  spac¬ 
ing,  dome  thickness,  density  of  the  dome  material,  and  material 
properties  of  the  dome  was  determined.  It  was  found  that  the 
most  important  parameter  in  determining  dome  effects  is  the 
quantity  p  rnh,  where  p  is  the  density  of  the  dome  material, 
cu  is  the  angular  frequency  of  the  transmitted  beam,  and  h  is 
the  thickness  of  the  dome  material.  The  quantity  o  u)h  is  the 
dominant  term  in  the  expression  for  the  modal  impedance  of  the 
dome , 


ih  E  i  h3E  .  ,  1  h  n2E2 

-ip  ha)  +  — * —  + - r—  [n  -  2n  +  1]  +  — « - n— ^ - * — —  , 

tub  12  tub  U)b^(o  u)  b  -  n  E  ) 

s  s 


( 


where  Zn  is  the  ratio  of  acoustic  pressure  to  normal  particle 
velocity  at  the  mean  surface  of  the  dome  for  the  nth  term  in 
the  sound  pressure  field  expansion.  E  is  the  elastic  modulus 

b 
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FIG.  2  -SOUND  PRESSURE  DISTRIBUTION  (DB)  FOR  A  TRANSDUCER 
RADIATING  INSIDE  A  CONCENTRIC  DOME  (NEARFIELD) 

ka  =  35.2  kb  =  43.9 

SHADING  =  COSINE 
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SOUND  PRESSURE  AMPLITUDE,  DB 


I 


FIG.  3  -SOUND  PRESSURE  DISTRIBUTION  (DB)  FOR  A  TRANSDUCER 
RADIATING  INSIDE  A  CONCENTRIC  DOME  (FARFIELD) 

Ka  =  35.2  kb  =43.9 

SHADING  =  COSINE 
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FIG  4  -  COMPARISON  OF  SOUND  PRESSURE  DISTRIBUTIONS  (DB)  FOR  A 
TRANSDUCER  RADIATING  WITH  AND  WITHOUT  A  CONCENTRIC 
DOME 
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FIG.  5 -COMPARISON  OF  SOUND  PRESSURE  DISTRIBUTIONS  (DB)  FOR  A 
TRANSDUCER  RADIATING  WITH  AND  WITHOUT  A  CONCENTRIC 
DOME 
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of  the  shell  ma terial ,  and  b  is  the  radius  of  the  mean  surface 
of  the  dome. 

Since  dome  effects  are  determined  principally  by  o  h, 

,  for  a  given  transmit  frequency  m,  it  is  possible  to  compact  the 

data  given  in  Ref.  2  to  sets  of  curves  which  demonstrate  the 

functional  dependence  of  the  performance  criteria  on  oh.  Some 

s 

typical  curves  are  shown  in  Figs.  6-10.  Figure  6  indicates  the 
dependence  of  the  two-dimensional  directivity  index  on  pgh. 

Figure  7  shows  the  variation  of  side  lobe  level  as 
a  function  of  pgh  and  Fig.  8  indicates  the  dependence  of  major 
lobe  width  on  p  h.  Results  for  dome  transmission  coefficient 

b 

and  dome  transfer  coefficient  are  shown  in  Figs.  9  and  10 
respectively.  The  transmission  coefficient  is  a  measure  of  the 
ratio  of  power  transmitted  through  the  dome  to  power  incident 
on  the  dome.  Dome  transfer  coefficient  is  a  measure  of  the 

» 

power  transferred  to  the  farfield  for  a  transducer  with  dome  as 
compared  to  a  bare  transducer.  This  latter  quantity,  a  measure 
>  of  the  acoustic  coupling  of  the  dome  and  transducer,  is  shown  as 

a  function  of  dome -transducer  spacing  (b-a)/\  due  to  its  strong 
dependence  on  the  spacing  parameter.  The  mathematical  definitions 

» 

of  the  performance  criteria  are  given  in  Ref.  2. 

C .  Viscoelastic  Domes 

o 

All  computed  dome  effects  in  the  studies  to  date 
have  assumed  that  the  dome  was  an  elastic  body.  The  application 
of  viscoelastic  materials  to  dome  construction,  for  example  the 
wire-reinforced  elastomeric  dome  for  the  AN/SOS-26  sonar  system, 
necessitate  the  investigation  of  the  effects  of  such  materials 
on  dome -transducer  interactions. 

rc - “ - 

It  is  possible  that  for  a  given  :u,  Zn  may  be  either  quite  large 

or  near  zero  for  a  fixed  n;  however,  these  conditions  are  rela¬ 
tively  unimportant  in  determining  the  farfield  of  the  dome- 
transducer.  Resonant  dome  conditions  are  discussed  in  detail 
in  Appendix  A. 
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FIG.  6  -  DIRECTIVITY  INDEX  VS  DOME  SURFACE  DENSITY  ft,h 


FIG.  7  —  SIDE 


FIG.  8  -  MAJOR  LOBE  WIDTH  VS  DOME  SURFACE  DENSITY  />sh 
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It  is  possible  to  modify  the  existing  dome -transducer 
model  to  include  viscoelastic  dome  materials.  This  modification 
requires  only  a  slight  change  in  the  form  of  the  differential 
equations  which  describe  the  motion  of  the  dome.  As  shown  in 
Appendix  B,  a  simple  viscoelastic  material  can  be  treated  in 
the  usual  manner  if  the  material  constants  are  complex  numbers. 
(The  material  constants  of  elastic  materials  are  real  numbers.) 

The  modal  impedance  of  a  viscoelastic  shell  is  similar  in  form 
to  Eq.  4.  However,  since  E  is  complex  for  the  viscoelastic 

O 

shell,  the  modal  impedance  has  both  7:eal  and  imaginary  parts. 
Losses  in  the  material  are  associated  with  the  real  part  of  Z  . 

The  modal  impedance  of  an  elastic  material  is,  of  course,  purely 
reactive  (imaginary)  as  shown  in  Eq.  4. 

The  only  practical  problem  associated  with  investi¬ 
gating  viscoelastic  dome  materials  is  the  unavailability  of 
realistic  values  for  the  complex  moduli.  The  'TRACOR  Rockville 
Laboratory  is  inves tigating  experimental  methods  to  determine 
these  moduli.  No  computations  of  viscoelastic  dome  effects  have 
been  executed  to  date  due  to  the  lack  of  realistic  material  con¬ 
stants.  However,  analysis  of  the  effects  of  a  rubber  dome  for 
the  AN/SQS-26  sonar  system  will  be  conducted  for  the  AN/SQS-26 
Project  Office,  Code  1631,  in  the  near  future.  The  rubber  dome 
and  typical  steel  domes  will  be  compared  in  terms  of  the  perfor¬ 
mance  criteria  previously  discussed,  using  the  best  estimates 
available  for  the  viscoelastic  moduli. 

D .  Nonconcentric  Dome  Shapes 

The  greatest  deficiency  of  the  present  dome -transducer 
model  is  the  unrealistic  dome  geometry.  The  next  logical  improve¬ 
ment  in  the  dome -transducer  model  is  the  replacement  of  the 
concentric  dome  with  an  elliptical  one.  It  is  contemplated  that 
the  analysis  will  proceed  as  in  the  past,  i.e.,  a  solution  to 
the  scalar  wave  equation  will  be  obtained  for  boundary  conditions 
expressed  on  the  transducer  and  dome  surfaces.  Since  these 


TRACOR,  INC. 


6500  TRACOR  LANE,  AUSTIN,  TEXAS  78721 


surfaces  are  expressed  in  different  sets  of  coordinates,  trans¬ 
formations  of  the  solution  forms  from  cylindrical  to  elliptical 
coordinates  will  be  required.  A  basic  problem  in  the  development 
of  these  transformations  is  the  computation  of  Mathieu  functions 
which  arise  in  the  general  solution  of  the  scalar  wave  equation 
expressed  in  elliptical  coordinates.  Methods  for  computing 
these  functions  generally  are  not  available. 

Possible  computational  methods  to  be  used  in  the 
generation  of  Mathieu  functions  have  been  investigated  during 
the  past  year,  and  successful  methods  have  been  developed  to 
the  point  that  the  solution  of  the  boundary  value  problem  appears 
feasible.  Further  work  on  obtaining  a  solution  to  the  bcundarj' 
value  problem  will  continue  in  the  studies  to  be  conducted 
during  the  following  year. 
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IV.  CONCLUSIONS 

The  studies  conducted  during  the  past  year  lead  to  the 
following  specific  conclusions: 

1.  The  dome  causes  extreme  nonuniformity  in  the  nearfield 
sound  pressure  distribution.  With  no  dome,  the  nearfield  sound 
pressure  distribution  is  relatively  smooth.  However,  the 
"nearfield  roughness"  as  defined  in  Ref.  2  increases  markedly 

if  the  transducer  is  enclosed  in  a  dome.  This  effect  can  increase 
the  radiation  loading  on  transducer  elements  as  well  as  decrease 
the  cavitation  limited  source  level  (through  the  development  of 
"hotspots"  in  the  nearfield) . 

2.  The  presence  of  a  dome  increases  the  side  lobe  level 
of  the  farfield  pattern.  Shading  of  the  transducer  elements  to 
obtain  reduced  side  lobe  levels  is  possible  with  a  dome;  however, 
the  side  lobe  level  reduction  appears  to  decrease  as  the  degree 
of  shading  increases.  Main  beamwidth  is  not  affected  signifi¬ 
cantly,  but  directivity  index  decreases  if  the  transducer  is 
enclosed  in  a  dome. 

3.  The  energy  transmitted  to  the  farfield  is  attenuated 
by  the  dome.  The  attenuation  is  due  primarily  to  an  interaction 
between  the  dome  and  transducer  and  not  to  losses  in  the  dome 
material.  Reflections  from  the  dome  change  the  radiation  loading 
at  the  transducer  face  and,  consequently,  can  affect  the  energy 
transferred  to  the  surrounding  medium.  This  interaction  is 
dependent  on  the  spacing  of  the  dome  and  transducer  and  it  may 

be  possible  to  reduce  the  interaction  by  proper  spacing  of  the 
dome  and  transducer,  at  least  for  limited  directions  of  transmission. 

4.  Degradations  of  typical  dome -transducer  performance 
criteria  such  as  directivity  index  depend  principally  on  the 
surface  density  (product  of  mass  density  and  thickness  of  the 
dome).  As  the  surface  density  increases,  degradations  in  the 
performance  criteria  increase.  The  surface  density  can  be  shown 
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to  be  the  dominant  factor  in  the  modal  impedance  of  the  concentric, 
unstiffened  dome.  It  is  possible  that  one  or  two  modes  of  the 
dome  may  be  in  a  resonant  or  antiresonant  condition  for  any 
given  transmit  frequency;  however,  these  modes  do  not  control 
the  effects  of  the  dome  on  transmitted  beams. 

The  above  conclusions,  again,  are  based  on  an  elastic, 
homogeneous  dome.  Other  studies  indicate  that  dome  stiffening 
members  can  further  degrade  the  acoustic  performance  of  domes." 

In  addition  to  studies  of  the  concentric  dome  model, 
an  extensive  investigation  of  methods  for  computing  Mathieu 
functions  has  been  completed.  These  functions  are  required  in 
order  to  treat  an  elliptical  dome  —  a  considerable  improvement 
in  dome  geometry.  Successful  methods  have  been  derived  and  the 
solution  of  a  boundary  value  problem  appears  feasible . 

Finally,  a  method  of  treating  viscoelastic  materials 
has  been  investigated.  A  simple  modification  of  the  differential 
equations  describing  shell  motions  is  necessary  in  order  to  include 
viscoelastic  effects.  The  modification  entails  the  replacement 
of  real  material  constants  (elastic  body)  with  complex  constants 
(viscoelastic  body) .  The  computation  of  viscoelastic  effects 
may  be  a  problem  due  to  the  scarcity  of  realistic  material  con¬ 
stants.  The  method  of  computation  is  simple,  however. 
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V.  RECOMMENDATIONS 

Continued  analytical  studies  should  be  directed  princi¬ 
pally  toward  improving  the  dome -transducer  model  in  terms  of 
more  realistic  geometry.  Additionally,  methods  for  treating 
stiffened  shell  structures  warrant  further  studies.  The  desi.r- 
ability  of  experimental  work  coordinated  with  analytical  work 
must  be  emphasized  in  an  exploratory  development  program  for 
sonar  dome  studies. 

A  three -year  program  of  analytical  and  experimental 
studies  to  improve  the  understanding  of  dome -transducer  inter- 
actions  has  been  submitted  to  BuShips ,  Code  1622E.  The  basic 
features  of  this  proposal  include  the  following: 

1.  Studies  to  improve  the  geometry  of  existing  dome- 
transducer  models, 

2.  Extension  of  studies  of  dome -transducer  interactions 
to  spherical  and  conformal/planar  configurations , 

3.  Related  experimental  studies  to  provide  verification 
of  analytical  results  and  assistance  in  generalizing  analytical 
models . 

In  particular,  studies  to  be  conducted  under  the  first 
year  of  the  program  include  the  development  of  a  model  consisting 
of  an  elliptical  dome  and  cylindrical  transducer  (for  which  some 
fundamental  research  is  completed) ,  studies  of  mathematical  tools 
which  have  application  to  the  interaction  problem  independent  of 
system  geometry,  and  experimental  measurements  of  the  vibration 
characteristics  of  stiffened  and  unstiffened  shells.  These 
experiments  will  suggest  methods  of  handling  stiffened  shells 
mathematically . 

These  studies  are  essential  to  a  fundamental  understand¬ 
ing  of  dome -transducer  interactions  and  the  subsequent  development 
of  improved  design  guidelines  for  sonar  domes. 


-a  — w ■.  waowyg/B  *  ~»--y  . 
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APPENDIX  A 


The  sound  pressure  field  outside  the  concentric  dome  is 
given  by  the  expression 


P2 


ft 

n 


(ka)  +  znanCJn(kb)  Nn'(ka)  -  Nn'(kb)  Jn'(ka)  ] 


H^(kr)  cos  n9  e  (2) 

where  the  symbols  are  defined  in  Sec.  II.  As  discussed  in  Sec. 

II,  if  Z  ,  the  modal  impedance  of  the  dome,  were  zero  for  all 

values  of  n,  the  dome  would  not  affect  transmitted  beams. 

However,  Zn  is  not  equal  to  zero  except  possibly  for  one  mode 

(one  value  of  n)  .  It  is  also  possible  that  Zn  -  00  for  one  value 

of  n,  in  which  case  the  amplitude  of  that  particular  mode  (term) 

becomes  zero.  These  two  distinct  features  of  Z„  and  the  resulting 

n 

effects  on  transmitted  beams  are  discussed  in  this  appendix. 

The  modal  impedance,  Zn ,  of  the  shell  is  given  by 
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where  u)  is  the  angular  frequency,  o  is  the  mass  density  of  the 

S  r» 

shell,  h  the  thickness  of  the  shell,  E_  ■  - E  is  Young's 

s  l-vz 

modulus,  v  is  Poisson's  ratio,  and  b  is  the  mean  radius  of  the 
shell . 


The  behavior  of  the  impedance  equation  exhibits  two 
distinct  features.  At  the  angular  frequency  given  by 


(5) 
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Zn  -*  «,  an  antiresonance  condition.  This  condition  results  since 
the  denominator  of  the  last  term  of  Eq.  4  reduces  to  zero.  At 
angular  frequencies  given  by 


Zn  =  0,  a  resonance  condition.  Equation  6  is  obtained  by  setting 
Zn  in  Eq.  4  equal  to  zero  and  solving  for  uuQ.  The  solution  of 
Eq.  6  has  two  real  roots,  hence  there  are  two  frequencies  for 
which  Zn  =  0. 


The  angular  frequency  associated  with  Z  -  differs  for 

^  Eq  2* 

each  mode,  (value  of  n)  but  has  a  unique  number,  ( — ^r)  ,  associ- 

°sb 

with  each  choice  of  the  parameters.  There  are  two  angular 


frequencies  associated  with  Z^  -»  0  for  each  choice  of  the  above 
parameters.  For  a  fixed  value  of  n,  the  angular  frequency  for 


which  Z  "*  °°,  lies  between  the  two  angular  frequencies  for  which 
Zn  =  0.  As  the  modal  number  increases,  both  zero -impedance 
angular  frequencies  increase,  but  they  are  never  equal  nor  is 


either  one  equal  to  the  angular  frequency  for  which  ZR  is 
undefined.  The  behavior  of  Z-^j  a  typical  modal  impedance,  is 
shown  as  a  function  of  ui  in  Fig.  11. 


The  question  arises  as  to  what  distortion  of  the  trans¬ 
mitted  beam  would  occur  if  the  operating  frequency  f  =  2mu  were 

chosen  close  to  ui  /2rr  or  oj  /2n  for  some  particular  mode  of  the 

o  00 

dome?  In  order  to  examine  this  situation  id  and  uj  were  deter- 

o  ® 

mined  for  n  =  12  and  typical  values  of  the  parameters  o  ,  h,  and  E  . 

S  b 
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The  value  of  tu  I  =  2.11694(10)^  rad/sec  for  which  Z-,~  -  *  was 
ro  n=12  12 

obtained  from  Eq .  5.  Two  values  of  u>  were  obtained  from  Eq .  6. 

/  O  A 

These  were  oj _ |  =  2.124474(10)  rad/sec  and  1.89007(10)  rad/sec, 
°  n=12 

for  which  Z ~  0.  Computations  of  farfield  patterns  using  the 
values  id  =  2.11694(10)^  rad/sec  and  uj  =  2.12447^(10)^  rad/sec 
were  carried  out  and  the  results  are  shown  in  Figs.  12  and  13. 

As  shown  in  the  figures  no  essential  difference  in  the  patterns 
is  observed.  (Any  differences  likely  are  due  to  the  difference 
in  ka  of  the  transducer  rather  than  modal  resonance  or  antireso¬ 
nance  effects.) 


In  addition  to  the  farfield  patterns,  modal  amplitudes 
were  obtained  for  0  £  n  £  100  for  the  two  frequencies  <uo  and 
tu^.  These  are  shown  in  Figs.  14  and  15.  In  Fig.  14  it  is  seen 
that  the  modal  amplitude  for  n  =  12  is  zero  as  expected,  since 

12 


Zin  -  m.  For  Z-^2  =  0  (Fig-  15),  the  modal  amplitude  reduces 


to 


12 


rfo  (ka) 


as  discussed  in  Sec.  II. 


12 


From  these  results  it  can  be  concluded  that  driving  the 
shell  at  a  resonant  or  antiresonant  frequency  for  any  one  mode 
does  not  distort  the  beam  unduly,  i.e.,  more  than  would  be 
expected  at  a  nearby  frequency.  This  result  is  due  to  the  fact 
that  only  one  or  two  modes  can  be  in  a  resonant  or  antiresonant 
condition  for  any  given  angular  frequency  tu.  (It  may  happen 
that  tura  for  a  low  order  mode  may  be  nearly  equal  to  the  lower 
value  of  tu  for  a  higher  order  mode.);  and  this  mode,  bounded 

O 

in  magnitude  between  zero  and  -y? - ,  does  not  control  the 

H“  (ka) 

total  response  of  the  dome. 
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FIG.  12- FAR  FI  ELD  SOUND  PRESSURE  AMPLITUDE  (DB)  VS.  RELATIVE 
BEARING  FOR  Z,2**c*=>(«= 2.1 1694  X  I04  RADIANS/SEC) 

u  =  0.3  E=3XI07  PSI  p  =7.355  X  |0~^— RCE'S^c2 

s  lN4 
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FIG.  13-FARFIELD  SOUND  PRESSURE  AMPLITUDE  (DB)  VS  RELATIVE 


BEARING  FOR  Z,~  =  0(^=2.124474  X  I04  RADIANS/SEC) 


u  =  0.3 


E  =  3  X  107  PSI 
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MODE  NUMBER  N 


FIG.  14-  NORMAL  MODE  AMPLITUDES  OF  SHELL  VS 

MODE  NUMBER  N  (NORMALIZED  TO  FUNDAMENTAL 
MODE)  Z|2  — »-oo 
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FIG.  15-NORMAL  MODE  AMPLITUDES  OF  SHELL  VS 

MODE  NUMBER  N  (NORMALIZED  TO  FUNDAMENTAL 
MODE)  Zl2  =0 
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APPENDIX  B 

This  appendix  presents  a  method  of  accounting  for  energy- 
dissipation  within  a  viscoelastic  dome  material .  The  equations 
of  motion  for  the  shell  model  of  the  dome  are  based  on  a  stress 
strain  relation  of  the  form, 

a  «  Es  e  ,  (7) 

where  e  is  the  longitudinal  strain,  a  is  the  longitudinal  stress, 

and  E  is  a  stiffness  parameter.  For  the  case  of  homogeneous, 
s 

isotropic,  elastic  shells,  E  is  called  the  plate  modulus  and 

s 

can  be  written  in  terms  of  the  Young's  modulus  E  and  the  Poisson 

o 

ratio  v  of  the  shell  material,  E  =  E/(l-v  ). 

o 

Energy  losses  in  sound  transmission  through  a  medium  are 
explained  in  terms  of  a  viscous  dissipation  mechanism.  This 
phenomenon  is  included  by  adding  one  term  to  the  stress  strain 
relation, 

a  =  Es  s  +  u  at  >  (?) 

where  u  is  the  viscosity  coefficient  and  t  is  time.  This  depen¬ 
dence  of  the  stress  on  the  strain  and  the  strain  rate  describes 
the  simplest  type  of  viscoelasticity.  For  steady  state  vibration 
with  angular  frequency  u),  the  stresses  and  strains  can  be  considered 
to  have  an  e~luJt:(i  =  /-T)  time  dependence.  Then  Eq .  8  takes  the 
form 

a  =  [Eg  -  iuHije  .  (9) 

Equation  7  and  Eq .  9  have  the  same  form.  The  correspondence  of 
Eq.  7  and  Eq .  9  suggests  that  the  shell  equations  for  a  viscoelas¬ 
tic  shell  may  be  obtained  from  those  for  an  elastic  shell  by  simply 
making  E  complex.  The  validity  of  using  a  complex  elastic  modulus 

b 

in  thin  plate  theory  is  discussed  in  Ref.  6. 

The  Fliigge  equations  for  a  cylindrical  shell  are  differen¬ 
tial  equations  for  the  shell  displacements  as  functions  of  the 
angle  9,  which  varies  from  -rr  to  rr  on  the  shell  surface. 
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Let  the  radial  and  tangential  displacements  from  equilibrium 

of  the  middle  surface  of  the  shell  be  denoted  by  u  and  u A , 

-j  r  fj 

respectively.  Fliigge's  shell  equations  are7 


where  h  is  the  thickness  of  the  shell  and  b  is  the  radius  of  its 
middle  surface,  p  is  the  mass  density  of  the  shell  material  and 
p  and  p+  are  the  pressures  on  the  inside  and  outside  shell 
surfaces,  respectively. 

These  equations  determine  the  motion  of  a  viscoelastic 
shell  with  the  complex  value  of  E  defined  by  Eq .  9.  The  deter- 

o 

mination  of  the  sound  pressure  field  produced  by  a  transducer 
radiating  inside  a  concentric  viscoelastic  shell  is  then  obtained 
in  a  manner  exactly  analogous  to  that  for  an  elastic  shell. 
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